A photonic cluster state machine gun 
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We present a method to convert certain single photon sources into devices capable of emitting large 
strings of photonic cluster state in a controlled and pulsed "on demand" manner. Such sources would 
greatly reduce the resources required to achieve linear optical quantum computation. Standard spin 
errors, such as dephasing, are shown to affect only 1 or 2 of the emitted photons at a time. This 
allows for the use of standard fault tolerance techniques, and shows that the photonic machine gun 
can be fired for arbitrarily long times. Using realistic parameters for current quantum dot sources, 
we conclude high entangled-photon emission rates are achievable, with Pauli-error rates per photon 
of less than 0.2%. For quantum dot sources the method has the added advantage of alleviating the 
problematic issues of obtaining identical photons from independent, non-identical quantum dots, 
and of exciton dephasing. 
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The primary challenge facing optical quantum compu- 
tation is that of building suitable photon sources. The 
majority of effort has been directed at single photon 
sources. Four single photons can be used in an inter- 
ferometer to produce a maximally entangled Bell pair 
of photons PP, and given a source of Bell pairs it is in 
principle possible to fuse them [2] into larger so-called 
cluster states |3j. These somewhat magical quantum 
states can be used for performing quantum computation 
via the simple procedure of making individual (single- 
qubit) measurements on the photons involved. Recently 
a promising new approach has been to produce Bell pairs 
directly [H H] via a radiative cascade in quantum dots. 
However, even an ideal such source would only reduce 
the overall resources required for a full optical quantum 
computation by a small factor. 

We will show that with current technology it is possible 
to manipulate certain single photon sources, in particular 
quantum dots, so as to generate a continuous stream of 
photons entangled in long strings of (various varieties of) 
1-dimensional cluster states. Using these strings, cluster 
states capable of running arbitrary quantum algorithms 
can be very efficiently generated by fusion. We analyze 
all error mechanisms and show that the error rates can be 
very low - close to fault tolerant thresholds for quantum 
computing - even if the source is operated for timescales 
much longer than the typical decoherence times. 

We begin with a highly idealized description of the 
proposal. Consider a source with a degenerate spin 1/2 
ground state manifold. The basis |f), |J.) denotes the spin 
projection along the z axis. Furthermore, imagine that 
optical transitions at frequency luq are possible only to 
a doubly degenerate excited state manifold. The excited 
states 1 1r) , | 4) have J z — ±3 /2H, thus only the (single 
photon) transitions |f) | f|~) and | J.) <-» | -IJ.) are allowed. 
Such transitions are well known to occur, for example, 
in quantum dots (QDs) which emit single photons via 



charged-exciton decay [B]. We only consider the emitted 
photons propagating along the z axis. Therefore, if the 
initial state of the source is ||) an excitation to 

the state | -ft) (| 1|)) followed by radiative decay, results 
in the emission of a single right (lcft)-circularly polarized 
photon \R) (\L)) and leaves the source in the state ||) 
Now, consider the initial state | f ) + \\), and a 
coherent excitation pulse with a linear polarization along 
the x direction. (The exciting pulse itself need not nec- 
essarily propagate along the z direction, which is useful 
for separation of the coherent and emitted light). Such 
a pulse couples equally to both transitions. Therefore, 
the processes described above happen in superposition, 
and the emitted photon will be entangled with the elec- 
tron: the joint state of both systems would be the Bell 
pair |T)i2) + Repeating such a procedure would 

produce GHZ-type entangled states, which are not use- 
ful for quantum computing, and for which disentangling 
the photons from the electron spin is difficult. Moreover, 
the GHZ state is highly vulnerable to decoherence. By 
contrast, the cluster states suffer none of these problems. 

To see how to create cluster states, we now imagine 
that before the second excitation of the system, when the 
state of the spin and the first photon is + 
the spin undergoes a 7r/2-rotation about the y-axis. Un- 
der this operation, described by exp(— iYn/A) the state 
evolves to (|T) + \Ri) + (- |T) + \Li) . A second pulse 
excitation, accompanied by a second photon emission, 
will now result in the two photons and the electron spin 
being in the state (|T)|-R 2 ) + \l)\L 2 ))\Ri) + HT)|i? 2 ) + |4 
)\L 2 ))\Li). In terms of abstract (logical) qubit encodings 
we will take \R) = |0), \L) = — 11). It can be readily 
verified that rotating the spin with another n/2 rota- 
tion, now leaves the spin and two photons in the state: 
|000) + 1001) + 1010) - |011) + 1100) H- 1 101) - |110) + 1 111) , 
which is exactly the 3 qubit linear cluster state. Repeat- 
ing the process of excitation followed by 7r/2 rotation, 
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FIG. 1: (Color online) A quantum circuit readily verified 
to output linear cluster state. For mapping to the clus- 
ter state machine gun, the top qubit line is the electron 
spin, the Hadamard gates are replaced by single qubit uni- 
taries exp(— inY/4) (requiring the careful tracking of certain 
phases), and the physical process of creating a photon with 
left/right circular polarization conditioned on the state of the 
electron spin becomes the controlled not gate which leaves 
the qubit (photon) in state |0) (i.e. \R)) if the electron spin is 
in state |0) (i.e. ||)), but otherwise flips it. Crucially, as de- 
picted, a Pauli Y error on the spin localizes; i.e it is equivalent 
to Y and Z errors on the next two photons produced. 



will produce a third photon such that the electron and 
three photons are in a 4-qubit linear cluster state. The 
procedure can, in principle, be repeated indefinitely, pro- 
ducing a continuous chain of photons in an entangled lin- 
ear cluster state. Note that one advantage of producing a 
cluster state is that the electron can be readily disentan- 
gled from the string of entangled photons, for example 
by making a computational (\R), \L)) basis measurement 
on the most recently created photon. In fact, since in 
general the initial state of the spin will be mixed, such a 
detection of a photon in state |i?)(|L)) polarization can 
also be used to project the spin to the | f ) (|J.)) state, 
and initializes the cluster state (either outcome is ok) . It 
can be readily verified that the whole idealized procedure 
just described is equivalent to the qubit quantum circuit 
depicted in Fig. [I] 

A general analysis of how cluster states are gener- 
ated by evolution of atoms in cavities undergoing general 
pumping and decay can be found in [7J, and interest- 
ing cavity QED proposals can be found in [5]. We will 
primarily focus on a specific implementation of our pro- 
posal, namely photon emission from a quantum dot, via 
the process of creation and subsequent decay of a charged 
exciton (trion). In practise the expressions we derive, 
such as the structure of the emitted photon wavepackets, 
can be easily applied to any systems which obey simi- 
lar selection rules, and the imperfections we discuss are, 
for the most part, generic. The importance of the selec- 
tion rules arises as follows. In semiconductor quantum 
dots, the J = 3/2, J z = ±1/2 states are naturally split 
off from the J = 3/2, J z = ±3/2 ones primarily due to 



confinement. They correspond to trions containing two 
electrons in the singlet state and a light hole or heavy 
hole respectively. We can consider only the heavy trions 
and neglect the mixing between them. In other systems, 
while the transitions to J z — ±1/2 may be energetically 
split off by an external field, or may simply have differ- 
ent couplings, gcncrically they will still lead to imper- 
fections equivalent to non-orthogonality of the emitted 
photons. Moreover processes in other systems tend to 
be slower, and temporally longer pulses may well also be 
required because of nearby energy levels. Although these 
problems can be remedied somewhat by applying proper 
filtration protocols to the output cluster state (at the ex- 
pense of larger loss rates) we focus on quantum dots for 
which the suppression is essentially perfect, the processes 
are fast, and the energy levels well separated. 

Although other options exist, we will consider from 
now on the situation where the 7r/2 rotations on the spin 
are performed by placing the quantum dot in a constant 
magnetic field of strength B which is directed along the 
y-direction (i.e. in the plane of the dot). The spin preces- 
sion at frequency lob — ge^B/h in the z-x plane therefore 
implements the desired rotation every T cyc i = tt/2ub- 
Suitably timed strobing of the dot by the excitation pulse, 
followed by the rapid exciton decay, will therefore enable 
the machine-gun-like generation of Id cluster state de- 
scribed above. 

The potential imperfections to be considered are as fol- 
lows: (i) The non-zero lifetime of the trion Td cca y means 
that the magnetic field causes precession of the electrons 
during the emission process. This leads to errors induced 
on the quantum circuit of Fig [I] however we shall find 
that they can be understood as implementing an error 
model on the final output cluster state which takes the 
form of Pauli errors occurring with some independent 
probability on pairs of (photonic) qubits. (ii) Interaction 
of the electron spin with its environment results in a non- 
unitary evolution of the spin. This evolution consists of 
two parts: decoherence (in which we include both de- 
phasing and spin flips) and spin relaxation. Decoherence 
is characterized by a Ti time. Fortunately we will see 
that both these processes also lead only to errors occur- 
ring independently on two (photonic) qubits at a time. 
Efficient cluster state quantum computation can proceed 
even if every qubit has a finite (though small) probability 
of undergoing some random error [§] . This implies that 
the protocol's running time is not limited by T2, while the 
errors are amenable to standard quantum error correction 
techniques for cluster states. Spin relaxation is charac- 
terized by a T\ time, and is a process which projects 
the spin to the ground state. In semiconductor quan- 
tum dots T\ times are extremely long T\ ^> T2 Td CC ay 
[TTj . Therefore, we shall not discuss the effects of this 
process further here. We point out, however, that it can 
be shown this process also leads to errors of a localized 
form, and so in principle is no obstacle to the continuous 
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operation of the device even for times much longer than 
T\ . (iii) The last source of error is related to the issue of 
ensuring the photons are emitted into well-controlled spa- 
tial modes. In practise this technological issue of mode 
matching (say by placing the dot in a microcavity) re- 
sults in some amount of photon loss error in the final 
state. Significant progress on this issue is being made for 
a variety of quantum dots [12j [13] , although we empha- 
size that for our proposal strong coupling to the cavity is 
not required. Fortunately photonic cluster state compu- 
tation can proceed even in the presence of very high (up 
to 50%) loss [TU], and we will not consider this source of 
error further. 

We now turn to detailed calculations of the error rate 
inflicted by imperfections (i) and (ii) discussed above. We 
first calculate the effect of a finite ratio of the trion decay 
time Tdecay time to the spin precession time. We denote 
by pn{j + tn) the state of the system (the quantum dot 
and photons) at time r after the n th excitation pulse, 
t n = nT cyc i e . By p n (t~) we mean the state of the system 
just before the n th excitation pulse (we assume the ex- 
citation is instantaneous). Following the excitation, the 
trion state decays, emitting a photon and leaving an elec- 
tron in the QD, the spin of which then precesses in the 
magnetic field. These lead to an evolution of the quan- 
tum state described by the following map (see [H] for 
details) : 

P (t n +r) = U\t){G + F)^p(t-)(G + F)U(r) (1) 

The unitary operator U = exp(iYu>BT) exp(iHoT) de- 
scribes the precession of the electron spin and the free 
propagation of the photons. The generalized creation 
operators & = C^ITXT I + G^U I, ^ f = 4U>U 
I — -^lITX! I; describe the excitation and decay process, 
adding a photon to the state. The trion states decay ex- 
ponentially with t, therefore we have omitted them from 
Eq. ([T]) (which describes the state of the system at times 
greater than the trion decay time, i.e. r» Td CC ay)- Note 
that the photons created in each cycle are well separated 
from the ones created in the previous cycles (formally, 
this is taken into account by the free propagation of the 
photons). 

Equation (JT|) describes a circuit isomorphic to the one 
in Fig. [l] The operator corresponds to a correct ap- 
plication of a CNOT gate. This happens with an am- 
plitude g{k), which depends on the photon's energy k: 

(ke\Gi\0)=g(k) = (fc ^gg^/4 - Here l**> = 4,J°> 
and e = L, R. The complex energy of the trion states 
is denoted by Z = ujq — iT/2, where Td CC ay = 1/r 
is their lifetime. The operator corresponds to a 
CNOT gate followed by a Y error on the spin qubit. 
This errored gate is applied with amplitude f(k), where 

(ke\Fm ^ fik) = i^ g %^ /4 . 

Let us for the moment treat the processes described 
by G^ and F^ as incoherent with each other. Then the 




FIG. 2: (Color online) Contour plots showing the probability 
of Pauli error on any given photon, as a function of g e fiB/T 
and (rr 2 ) _1 . A stronger field causes faster precession which 
increases a chance of the error during decay, but reduces the 
standard dephasing error due to a finite T%. The inset is a 
plot of the good mode function \g(k)\ (blue) and error mode 
function |/(fe)| (red) at g e jj,B/T — 0.15, from which it can 
be deduced that spectral filtering can reduce the error rates 
further. 

resulting state is described by the circuit of Fig. [I] with a 
probability p B = ||/|| 2 = 2(g %%f+ 2r2 , that each CNOT 
gate is followed by Y error on the spin qubit. As noted in 
Fig. 1, a state with a Y error on the spin after generation 
of the n th photon (i.e, after the n th CNOT), is equivalent 
to a state Y and a Z error on the n th + 1 and n th + 2 pho- 
tons, with no error on the spin. Note that the error prob- 
ability increases with magnetic field strength, because 
the spin can precess more during the lifetime of the trion 
[TS] . Therefore it is advantageous to consider relatively 
low magnetic fields, for which g e pB <C T. Taking the co- 
herence between G' and F^ into account, it can be seen 
that a unitary correction e lY< ^ with tan (j> = \ (g\f) / (g\g)\ 
yields a further improvement of the error rate. We also 
point out that as g{k) is more localized around ojq then 
f(k) (inset of Fig. selection of photons with energy 
| k — luq I > A would yield a lower error rate at the expense 
of (heralded) loss. 

The calculation above ignores the possibility of the ex- 
citon dephasing |16] during the decay process. Pure de- 
phasing, in which both excited levels evolve the same 
(random) phase, will have no affect on the entanglement 
in polarization with which we are concerned. Cross de- 
phasing (experimentally seen to be very small |16j ) will 
lead to Z-errors on the qubits, which also localize (See 
[2] for a detailed discussion). 

We now turn to the issue of decoherence of the spin as 
a result of its interaction with the nuclei in the quantum 
dot. Assuming Markovian dynamics (discussed further 
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in [H|), it is well known [T7] that the resulting dephas- 
ing and spin flip dynamics are equivalent to the action of 
random Pauli operations X, Y, Z with some probabilities 
Px,Py,Pz- The probabilities p x , p z are suppressed due to 
the presence of the magnetic field, while p y is charac- 
terized by T2, the dephasing time. This can readily be 
shown to give p y — i(l — e~ Tcycl °/ T2 ) as the probability 
of a given spin error in 1 cycle. We already noted that 
a Y error on the spin becomes a pauli error on the next 
2 photons. Similarly a Z error at the end of the n th cy- 
cle is equivalent to a Z error on the n th + 1 photon, as 
can be seen from Fig. [l] and the fact that the operator 

C l NOT(-Z'spin ® -fphoton) and (Igpin ® ^photon)CNOT have 

similar actions on the states 1 00) and |10). As X = iZY, 
an X error again affects only the next two photons to be 
generated. 

In Fig. [2] we plot the total probability of error on any 
given qubit, 1 — (1 — ps)(l — p y ), as a function of the 
two dimensionless parameters g e /iB/T and (rT 2 ) _1 . We 
include the aforementioned easily implemented unitary 
correction. Although free induction decay T2 may be rel- 
atively short in low magnetic fields, using a spin echo 
pulse at half way along the cycle time can extend Ti con- 
siderably, and remove the dephasing caused by a wide 
distribution of nuclear (Overhauser) magnetic fields (of- 
ten termed inhomogeneous broadening and characterized 
by a T| time). To estimate an achievable error rate, we 
consider a decay time of 1/T = 100 ps, and a dephas- 
ing time of Ti = 1 /is with the addition of the spin echo 
pulses (a lower bound of 3 [is have been measured in 
high magnetic fields [13]). This gives (FT?)- 1 = 10~ 4 . 
From Fig. [2j one can deduce that a probability of error 
less then 0.2% can be achieved by applying a magnetic 
field of 15mT (we take g e — 0.5). We note that even 
without the spin echo pulses, error rates of about 1% are 
achievable, which enables the production of considerable 
longer and higher quality optical cluster states than those 
produced by current methods. 

So far we have considered pulse excitations that are 
timed to coincide with (integer multiples of) 7r/2 rota- 
tions of the spin. In fact it can be advantageous to some- 
times wait for a full tt rotation to occur. This has the 
effect of emitting subsequent photons which are redun- 
dantly encoded [2J. Fusing together such qubits gives a 
highly efficient method for producing higher-dimensional 
cluster states which are universal for quantum comput- 
ing. Photons which undergo fusion can be spectrally fil- 
tered (via a suitable prism), such that if they fail to pass 
the filter they can still be measured and removed from 
the cluster state. This filtering does not lead to an in- 
crease in loss error rates, but simply decreases the overall 
success probability of the fusion gates. 

Current experiments produce photonic cluster states 
via spontaneous parametric downconversion [19 , and 
would seem to be limited to producing 6 to 8 photon clus- 
ter states. Our proposal in principle can produce strings 



of thousands of photons; however initial experiments will 
be limited by collection efficiency. With the parame- 
ters above a simple analysis shows that we would need 
a collection+photodetection efficiency of about 18% for 
a demonstration of on-demand 12-photon cluster states, 
where the full 12 qubits are expected to be detected about 
once every 10 seconds. 

Finally, our proposal is suggestive of an efficient mech- 
anism for entangling matter qubits [5J HH], and we feel 
this is a topic worthy of further investigation. 
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SUPPLEMENTARY MATERIAL FOR THE PHOTONIC CLUSTER STATE MACHINE GUN 

Note: In the published version this document accompanies the above paper as supplementary material in the form 
of an EPAPS archive - i.e it is the same material as can be found at reference [15] above. The references in this section 
are self contained and can be found at the end of this document. 



CALCULATION OF EMITTED PHOTON WAVEPACKETS 



In this section we compute the emitted photon wavepacket, taking into account the fact that the electron spin is 
precessing during the emission. This effect and the finite T 2 time of the electron, are the dominant source of imperfect 
cluster state production in the machine gun. Any notation not defined here is defined in the paper. 

The Hamiltonian we consider is 

H = Hq + Hi n t 

where H n = ui a P + Hb + Hem- Here P is a projector on the dot excited state (trion) manifold and Hem is the free 
EM field Hamiltonian. 

The Zccman interaction for the electron is Hb = g e fJ-BS e • B, where /i is the Bohr magneton, and g e is the effective 
gyromagnetic ratio of the electron in the QD. For the low magnetic fields we are dealing with, we can assume gh = 
for the heavy hole [1]. We abbreviate \g e \J^B — * 9e- We take B = By. The interaction Hamiltonian of the QD with 
the photon field, Hi nt is given in the rotating wave and dipole approximation by 



Hint = 5> fc (|T) <ft| 4, fe + 14) m a{ ik + h.c) 



where Vk are the coupling constants which depend on the details of the specific quantum dot in question, and the 
(heavy) trion states are denoted by | ft) , | JJ.) respectively. 
We want to calculate matrix elements of the form: 

(T, Rk\ Uit) \1t) = ^-.J dEeM-iEt) (T R\ Q{E) |ft) (2) 
Here \R k ) and \L k ) are left and right circularly polarised single photon states created by a) R k ,a L k , and 

is the resolvent of the hamiltonian. 

Denote by Q the projector on ground state manifold. Then the matrix block inversion formula gives 

Q5(E)P= = Q „ v 



E - QHQ E - PH P - PR(E)P 
with R(E) = Vt^qV. 

The operator R(E) can only couple the state |f|~) (or |JJ.)) to itself (recalling that the hole is assumed to not precess). 
Therefore 



E - PH P + PR(E)PJ , E - Lo - r{E) 



where r(E) = J dk E _ k , and s,s' = |ft) , |JJ.) . When performing the contour integral (1), we can replace r(E) with 
the lamb shift A and the inverse decay rate T. In the following we shall take Vk = ^T/2n for simplicity and absorb 
A into the definition of loq. This gives 

Q h (E) = 5 SS , t 



G 



Next, we consider 



Q e {E) = 



Q 



E - QHQ J a , a 
1 ( E-hk -ig e B 



E - hk ig e B 
-ig e B E - hk J , 



D e (E) \ ig e B E-hk ) a , a 
where D e (E) = (E - hk) 2 - g 2 e B 2 . We now have the final result: 

(a,e k \g(E)\s) = V k e G(E) h aa , =£ g(E)s'= 
Writing the matrix elements explicitly: 

(1,R k \g{E)\1t) = v k 

(l,R k \g(E)\it) = v k 

{\,L k \g{E)\$) = v k 

(l,L k \g{E)\H) = \), 



(E - hk) 

De(E)(E-cu + iT/2) 

igeB 

D e (E)(E-cu Q + iT/2) 

-ig e B 

'D^E)(E-cu Q + iT/2) 

(E - hk) 

: D^(E)(E-uj + iT/2) 



-(i,R k \g(E)\n) 
(\,R k \g{E) no 



We now consider the contour integral. The two poles contributing in the limit t — > oo are E = k—g e B and E = k+g e B. 
We denote Z = lu — iT/2 and a standard complex integration gives 

(t,i? fe | U(t) = [exp(-i fle Bt)/i +exp(ig e Bt)f 2 ] 



where /i = ( k +g^B-z )i h = ( k -gB-z )- Above and in the following we omit the phase factors e tkt . Using the 
notation 



g(k) = h(k) + f 2 (k) 
f(k) = .h(k)-f 2 (k) 



v/r/27r(fc - Z) 
(k - Z) 2 - (g e B) 2 

V^J^geB 



this simplifies to 



Similarly 



which yields 



(k - Z) 2 - (g e B) 2 
C\,R k \ U(t) \t) = cos(g e Bt)g(k) - i sm(g e B)f(k). 
{i,R k \ U(t) |fr) = iexp(-ig e Bt)f! - i exp(ig e Bt) f 2 



(i,R k \ U(t) |ft) = i [-ism(g e Bt)g(k)+cos(g e Bt)f(k)} 
= sm(g e Bt)g(k) + icos(g e Bt)f(k) 

An analogous calculation shows that the amplitudes for starting with the down trion are 

(T,i fc | U(t) \H) = - (i,R k \ U(t) \ii) = -sm(g e Bt)g(k) - i cos(g e Bt)f(k) 

and 

(l,L k \ U(t) |^) = (T,i? fe | U(t) |ft) = (tR k \ U(t) |ft) = cos(g e Bt) 9 (k) -ism(9eB)f(k) 

Starting from the ground state manifold, the excitation and subsequent photon emission and spin precesion can be 
broken into "good" and "bad" parts, by introducing the operators: 
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G f = G^|T))(T|+G t L ||)a|, 

^ = 4u>ui-*iit>ui> 

The operator corresponds to a correct application of a CNOT gate. This happens with an amplitude g(k), 
which depends on the photon's energy k: (fce|G||0)) = g(k), the operator corresponds to an errored gate which is 
applied with amplitude f(k), where (ke\F^\0) = f(k). 

Starting from the state | f), we get 

\ip) = U^ Rk , t \ tRk) +U URktt \ l,R k ) 
which we break \ip) into the good and bad states: \ip) = \ip g0 od) + \ipbad) as follows: 



(T,-Rfc| "4>good) = cos (g e Bt) g(k) 
(l,Rk\ipgood) = sin(g e Bt)g(k) 

and 



—i sin (g e Bt) f(k) 
zcos (g e Bt) f(k) 



(\,R k \ Vw) 
(l,Rk\ Vw) 

Note that \ip bad ) = j^Y spm \ip good ) where Y spm ||) 
and F is 



i \ l) and Y spin |J.) = — i ||) . The relation to the operators G 



(hRk\i> 3 ood) = (hR k \u(t)G\V 

(LRkWgood) = {l,Rk\U(t)G\1) 

and 



(t,i4|Vw) - (T,i? fe |t/(i)F|t> 

(l,i? fe |^w) = (1,R k \U(t)F\1) 

where U(t) denotes the free propagation. 

Similarily, if we start from the state | J.), the resulting state would be 

M = Ui, Lk M ],L k )+U LLkA \ [,L k ) 

Lets break into good and bad states: = \ip goo d) + \ipbad) as follows: 

(T,£fc| ipgood) = - sin (g e Bt) g(k) 
(i,L k \ ipgood) = cos (g e Bt) g(k) 

and 

(T,£ fe |Vw) = -i cos {g e Bt) f(k) 
(l,£fc|Vw) = -i sin (g e Bt) f{k) 

The relation to the operator G and F is similar to the one described above. 
Note that again |Vw) = j^rfeY spin \ip good ) ■ 
If we start from the state | |) + | I), the resulting state would be 



|V> = Cv(ffeBt)(5(A:) + /(fc)y ap i„)(| T,«fc> + U,£fc» 
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where U y denotes spin rotation around the y axis. The bad part of the wavefunction is therefore equivalent to a 
Pauli Y error on the ideal cluster state - and as discussed in the paper this error on the spin can be turned into errors 
on the either the two photons emitted the error occurs on the spin, or as Z errors on the photon which has already 
been emitted and the first emitted photon. This latter possibility shows that this "quantum scapegoat" effect can act 
backwards in time! Note that while the error is coherent, the coherence between \ip g0 od) and \ipbad) can be removed 
by applying stabilizers randomly to the photons. In fact the good and bad wavepackets have non-trivial overlap, and 
in the next section we consider how our knowledge of this can allow us to apply a simple unitary correction which 
increases the amplitude of ideal cluster state produced. 

UNITARY CORRECTION 

Considering the "bad" amplitude f(k) we can clearly write it as 

f(k) = (f(k)-ag(k)) + ag(k) 
with a = (g\f)/(g\g)- Note that first term is orthogonal to g(fc).Then 



H>) = ^=U y (g e Bt)[(l + aY spin )g(k) (| T,i? fc ) + | |,L fe )) + (f(k) - ag(k))Y spm (| T.ifc) + | i,L k ))] 
(Note that a is purely imaginary.) 

We can now make the correction exp(— i(f>Y sp i n ) where tan</) = a. Figure 2 of the paper is plotted assuming this 
simple form of unitary correction has been carried out. 



EXCITON DEPHASING 



In the paper we only briefly mentioned the spectral dephasing which will occur while the system is excited. Our 
intuition contrasted with that of others, namely we felt that this process would not affect the entanglement of the 
state - in particular with respect to the polarization degrees of freedom we are interested in - but would only lead to 
the emitted photon wavepackets being in a mixture of different frequencies. This in turn would only affect the (small 
fraction) of photons which have to go through fusion gates, and such photons can be filtered before entering the gates 
in a way which will only lead to a change in the success probability of the (non-deterministic) gate. That is, such 
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filtering need not even lead to a loss error (as explained below). As such the only effect will be that we need to use 
more photons - but the overhead is some constant factor. 

Nevertheless, to be sure of how the device behaves and of the potential magnitude of this effect, we turn now 
to the necessary calculations. We begin with a description of the results and a heuristic discussion. There are two 
excited state (exciton) energy levels to consider, which in the paper we denote |^). These levels are degenerate. 
We consider both pure dephasing and cross dephasing. In pure dephasing of these levels, their couplings to the 
environment are identical - in practise this means that both states could evolve the same (random) phase. Note 
that this will not affect the relative phase between the two states (such a relative phase would ultimately become a 
relative phase between the two terms in the entangled photonic cluster state). Therefore, it will have no bearing on 
the quality of the entanglement in the polarization degrees of freedom. It will, however, cause the emitted photons to 
have a broader range of frequencies (broadened linewidth). In [2] such pure dephasing was measured, and from their 
results one can infer that the new linewidth would be approximately double the natural linewidth. However, we stress 
that this broadening is not of particular importance to the operation of our protocol, for two reasons. Firstly the 
majority of the emitted photons will never go through any optical element where (for instance) they may be required 
to interfere with other photons. Rather they will simply be measured directly, and this polarization measurement 
will not care about their frequency. For the photons that do need to undergo Type-II fusion gates (to cross link with 
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other qubit lines) the incoherence over the emitted frequencies can be removed (if necessary) by suitable spectral 
filtering. This filtering can be done in an "active" manner with a suitable prism wherein any photon which does not 
pass through the filter is simply diverted into another path and measured in an appropriate basis - in terms of the 
cluster state this simply removes the qubit. As such the only effect of such imperfections is to change the probability 
of success of the Type-II fusion gate. However efficient quantum computation is obviously possible regardless how 
small this probability is, as long as it is finite. 

Let us now show some calculations to bolster the above discussion. We treat the pure exciton dephasing as a 
Markovian process and describe it using a Lindblad operator L = \/7<i/2(|ft) (ft | + | JJ-)(JJ- 1) where we use the same 
notations as in the paper. The density matrix of the system evolves according to a standard master equation. We 
compare the probabilities described by this density matrix to the probabilities |(/(fc)| 2 corresponding to no error, and 
\f(k) |, corresponding to a pauli error, which were calculated in the paper. In the presence of pure dephasing, the 
probabilities corresponding to no pauli error are 

|.gdc P hasc(fc)| 2 = ^ J dtij ^ 2 e(- r / 2 -^/ 2 )* 1 e(- r / 2+ ^/ 2 )* 2 cos(6i 1 )cos(6i 2 ) e ife ^-* 1 ) 

+ h r dt * 1^ dt ie (- r / 2 -W2)t 2e (-r/2+ 7 ,/2) tl CQs{bti) cos^e^-t!) 

71 Jo Jo 

In the above b — g\iB is the zeeman splitting, and we have shifted k — > k — ui a . The off-diagonal matrix elements 
corresponding to pauli errors have the same form as Eq. but with sin functions replacing the cos functions. Although 
the function |<7dephase(&)| 2 is wider then \g(k)\ 2 , the total probability is equal to the total probability without dephasing 

J ' dk\g dcphasc {k)\ 2 = j dk\g(k)\ 2 , (3) 

for any value -fd- 

In Figure 1 we plot the probabilities |gdc P hase(fc)| 2 versus \g(k)\ 2 . The parameters chosen for the plot are b/T = 0.15 
(as in the inset of Fig. 2 of our manuscript) and 7d/r = 1 [1], which of course enters only into |gdcphasc(&) | 2 - Likewise, 
Figure 2 plots |/ dc pha S c(fc)| 2 vs. \f(k)\ 2 . In Figure 3, we compare \gde P h asc (k)\ 2 and |/d op haso(fc)| 2 showing that spectral 
filtering is still possible also in the presence of pure exciton dephasing. 

A potentially more serious source of error would be cross dephasing. Experimentally such dephasing has not been 
observed, despite an attempt to measure it in [2], where they could only lower bound the dephasing to be at least 20 
times smaller than the pure dephasing. Even if this cross dephasing were exactly equal to the lower bound obtained 
in [2], the effect on the emitted photons would ultimately lead to Z-errors on the qubits, which also localize. This 
potential increase in the pauli error rate is negligible compared to the other sources of pauli error we considered in 
the paper. 

NOTE ON SPIN-BATH DYNAMICS 

To calculate the error probablilty on the photonic qubits we have assumed Markovian dynamics for the interaction 
of the spin with its environment. Previous studies [3-7] suggest a non-Markovian decay in time of the coherence 
in the reduced density matrix of a spin coupled to a nuclear spin bath. We note that for the short times we 
are interested in, the Markovian assumption actually overestimates the single qubit error rate. The resulting non 
markovian correlations of the errors on the different qubits will be studied in a forthcoming publication. However, we 
note that these correlations can be dealt with, see for example the discussion in [4]. 
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